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ABSTRACT 

Accretion rates onto AGN are likely to be extremely variable on short timescales; much 
shorter than the typical cooling time of X-ray emitting gas in elliptical galaxies and 
galaxy clusters. Using the Langevin approach it is shown that, for a simple feedback 
system, this can induce variability in the AGN power output that is of much larger am- 
plitude, and persists for longer timescales, than the initial fluctuations. An implication 
of this is that rich galaxy clusters are expected to show the largest and longest-lived 
fluctuations. Stochastic variations in the accretion rate also mean that the AGN in- 
jects energy across a wide range of timescales. This allows the AGN to maintain a 
much closer balance with its surroundings than if it was periodically activated. The 
possible non-linea r correlation between Bondi accretion rate and jet power, found by 
lAllen et all (|2006l ). can be explained if the instantaneous accretion rate, scaled by 
jet power, varies log-normally. This explanation also implies that the duty cycle of 
AGN activity i ncreases wi t h the radiative losses of the surroundings, in qualitative 
agreement with lBest et al.l ()2005l ). 
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1 INTRODUCTION 

It is well known that elliptical galaxi es are commonly th e 
hosts of powerful radio AGNs (e.g. iMcLure et all |2004 ). 
These sources give rise to lobes of radio emission embed- 
ded in the X-ray emitting gaseous haloes surrounding the 
host galaxies. There is also growing theoretical evidence 
that outflows from AGNs play a si gnificant role in the evo- 
lution of the i r surr o undings (e.g. Briiggcn & Kaiser 2002: 
Proton et all l200rj : ISiiacki fc Springeli l200rj : iBower et al l 



2006). 



In theoretical work, outflows from AGNs are com- 
monly used to prevent catastrophic radiative cool- 
ing i n both elliptica l gala xies (e.g. iTabor fc Binnevl 



19931: Binnev fc Tabor] 19951) and galaxy clu st ers (e.g 



Churazov et al.l 2001 : Briiggen fc Kaisei 20021: Briiggenl 



20031: iBasson fc Alexander! 120031: lOmma et all 



Dalla Vecchia et al. 20041 : ICattaneo fc Tevssierl 



However, for there to be a long term balance between the 
heating and cooling processes there must be some sort of 
feedback mechanism through which the cooling gas triggers 
AGN outbursts. The level of jet activity appears t o be a 
function of the environment and b lack hole mass (e.g. lBurnsl 
ll99d : lBest et al.ll2005l : [BestJl2007h . In addition, observations 
also suggest that the heating must be relatively gentle and 
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well matched to its environment (e.g. iFabian fc Sanders! 
l2006h . 

In numerical simulations, the jet power of an AGN is 
often related directly to the effective inflow rate of material 
through an inner boundary of the computational grid (e.g. 
IVernaleo fc Reynolds! 120061 : ICattaneo fc Tevssierl |2006| ). In 
reality the picture is not so simple since the overall accretion 
rate onto a supermassive black hole, located in an ellitical 
galaxy, is probably a c ombination of its local, and large- 
scale, env ironments (e.g. Sokcr 2006; Hard castle et al. I l2007l : 
ICenll2007l ). The local accretion rate is p robably governed by 
Bondi accretion (e.g. I Allen et al.l 1200a ). However, numeri- 
cal simulations suggest that Bo ndi accretion is not steady, 
but can be highly variable (e.g. Edgar 2004, and references 
therein). We also observe extremely powerful , cluster-scale 
outbursts, e.g. Hydra-A l|Nulsen et al.l 120051) . These out- 
bursts are probably triggered by accretion of material from 
the wider cluster environment, rather than Bondi accretion. 
In this scenario, one could envisage cold 'blobs' of material 
which drop out of the flow around a galaxy, 'diffusing' to- 
wards the central black hole. Clearly this process would be 
expected to be intrinsically variable. 

Most of the infalling material will have appreciable an- 
gular momentum; somehow this gas must be funnelled down 
onto the black hole, presumably through an accretion disc. 
This process will mediate the accretion rate from the exter- 
nal system that surrounds the black hole. Depending on the 
disc processes, the instantaneous black hole accretion rate 
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need not be strongly correlated with the external (systemic) 
accretion rate. 

Another source of accretion rate variability would be 
tidal interactions a nd mergers with nearby galaxies (e.g. 
Hopk ins et al]|2006h . Furthermore, the dynamic state of the 
accreted medium will also affect the black hole accretion 
rate. In this respect, some advances in numerical models 
have been m ade using initial cond itio ns based on t he sim u- 
lation of the ISpringel et all l|200ll ) bv lHeinz et al. (2006). 

In short, it is easy to understand how the accretion rate 
onto an AGN could be highly variable, on short timescales, 
compared to the approach assumed in many hydrodynamic 
simulations. However, adequately resolving the required 
range of scales is clearly impractical with current computing 
constraints: most numerical simulations of AGN feedback do 
not even resolve scales as small as the Bondi radius, let alone 
the accretion disc. 

Although there is not extensive literature, the impor- 
tance of random, or stochastic, accretion is slo wly becom- 
ing more apparent. Hopk ins fc Hernquisj (|2006l ) employ the 
phenomenon to investigate the fueling of low-level AGN, and 
its effect on the black hole - bulge relations. In their model, 
they assume that cold gas stochastically accretes onto a cen- 
tral supermassive black hole at a rate set by the dynamics 
of that gas. In particular, clouds of molecular gas collide 
with the supermassive black hole. This provides a distinc- 
tion between local, low-luminosity quiescent AGN activity, 
and violent, merger-driven brig ht quasars. 

iKing fc Pringlel i|2006h and lKing fc Pringlel (|2007l ) also 
study the effect of small-scale, randomly oriented accre- 
tion event, as a me chanism for feeding nearby AGN. 
IKing fc Pringlel (|2006l ) conclude that supermassive black 
holes can grow rapidly if most of the mass comes from a 
sequence of randomly oriented episodes whose angular mo- 
menta are no larger than the angular momentum of the black 
hole. This means that the black hole has a low spin, and 
therefore a low radiative efficiency, so rapid growth is pos- 
sible without exceeding the Eddington limit. Their model 
requires the accretion hydrodynamics to be chaotic at the 
lev el of the size of the disc: < O.lpc. This model is extended 
bv lNavakshin fc Kind (|2007T) who find that, due to the ran- 
dom orientation of the accretion disc in these events, energy 
can be fedback into the surroundings in an almost isotropic 
manner. This is important since it is much more efficient 
t ha n anisotropic heat i ng. 

iNipoti fc Binnevl (|2005T ) also argue that the mechanical 
luminosity of an AGN is expected to fluctuate across a wide 
range of timescales, while the X-ray luminosity of the cluster 
gas will only vary slowly and weakly with time. Since energy 
injection by AGN is expected to balance the radiative losses 
of the cluster gas, it makes sense to construct a dimension- 
less variable where the AGN luminosity, H, is scaled by the 
radiative losses, Lx- This eliminates variations in H which 
reflect the size of Lx - The authors argue that the probability 
density function of ln(H/Lx), at any time, is Gaussian; in 
other words, H/Lx is lognormally distributed, which is an 
example of a skewed probability distribution. 

Throughout physics, stochastic variati ons are often 
addressed using Lan gevin dynamics (e.g. ISchussI Il980l : 
iKarlin fc Ta^iorfll975l ). The Langevin equation is a stochas- 
tic differential equation in which force terms are added to 
Newton's second law to approximate the effects of neglected 



degrees of freedom. One term represents a frictional force, 
the other a random, fluctuating force. In this method, a ran- 
dom variable / represents physical processes that happen 
rapidly and which may be on smaller spatial scales than 
are resolved in the model. The principle assumptions for 
this fluctuating force are: / is statistically independent of 
the physical variable of interest, X. The variations of / are 
much more rapid than those in X, and the mean of / is zero. 
Therefore, / represents a source of 'white noise'. This prin- 
ciple can be applied directly to simulations of AGN feedback 
in elliptical galaxies and galaxy clusters. 

In this letter, we describe the use of the Langevin equa- 
tion in predicting the variability of AGN power output, and 
its dependence on environment. A more complex feedback 
model is used to explain the possible origin of the non-linear 
correlation between the Bondi accretion rate a nd jet power 
in a s ample of elliptical galaxies reported in lAllen et al.l 
(2006). Finally, we summarise the main findings of this work. 



2 GALAXY AND GALAXY CLUSTER 
VARIABILITY MODELS 

Observations suggest that elliptical galaxies and galaxy 
clusters must be rela tively close to steady state (e.g. 
iFabian fc Sandersfe oOG) with the mass of cold gas and new 
star formation lower than expected (|Edgell200li ). although 
see iRaffertv et all |2006). In addition the spectra also indi- 
cate the mass dropout rates are much lower than expected 
from classical estimates (e.g. IVoigt fc F abian 20(3). 

We can construct a very simple, but instructive, phe- 
nomenological model to investigate the effect of AGN feed- 
back. Let us assume that the accretion rate increases if the 
radiative losses exceed the energy injection rate, and vice 



H), 



(1) 



where H is the energy injection rate by the AGN, Lx are the 
radiative losses of the cluster atmosphere, and K\ is a con- 
stant. If we also assume that H and m vary on much shorter 
timescales than Lx, we can treat Lx as approximately con- 
stant. In the absence of stochastic accretion we can write 
H — K^rn, where K2 = r}(? ~ 10 20 cm 2 s -2 . Solving equa- 
tion JTJ), assuming that the accretion rate is zero at t = 0, 
gives, 



m(t) = — 



exp 



(2) 



where r = l/(KiK2)- In other words, the system tends to- 
wards a state with a steady accretion rate of Lx/ K2, as time 
passes. 

If there is a time-delay in the AGN response we would 
expect the accretion rate to overshoot the steady-state value, 
and subsequently perform damped sinusoidal oscillations 
around this point. Of course, if the time-delay is so large 
that the overshoot is equal to the steady-state accretion rate, 
then the system will be unstable. 

The value of K\ is slightly harder to evaluate; we esti- 



Rapidly varying accretion and A GN feedback 3 



mate it in the following way, 
dm rh 

d£ ^cool 



KiLx. 



(3) 



We can write Lx = 5kbTrh c / (2/j,m p ) , where m c is the mass 
dropout rate from the ICM, which is distinct from rh - the 
accretion rate onto the black hole. Thus we have 



T 



We m ight reasonably expect rh c /rh ~ 100 (e.g. IPope et al.l 
2006) so that r ~ 10 5 yrs. This is much shorter than typical 
radiative cooling times, thus justifying our initial approxi- 
mation that the accretion rate varied on comparatively short 
timescales. This also indicates that the heating response 
quickly establishes an equilibrium, but that r increases with 
the temperature of the ICM, suggesting that more massive 
clusters have a longer response time. 

Now, if we also consider the effect of stochastic accre- 
tion, we can write H — Kim + AH, where AH is a random 
variable that can be positive or negative. Since the system 
also has a steady-state accretion rate, (rh) = Lx/K2, we can 
write rh = (rh) + Am, 



tc 



dAm Am A „. . 

— =-—- Kl AH®. 



(•») 



Equation (JSJ) has the form of the Langevin equation where 
—KiAH — awf(t). fit) is a Gaussian 'white-noise' process, 
with zero mean and unit standard deviation; <tn, is the stan- 
dard deviation of the noisy process. The variance of Am is 
given by, 



(6) 



which demonstrates that seed fluctuations of (rh) j\Jr <SC 
{rh} can produce AGN power fluctuations that are compa- 
rable with the mean accretion rate. 

The transformation of a random variable through a dif- 
ferential equation, to ob tain a probability dist ribution, is 
notoriously complex (see IChatteriee et al. 120021 , for an ex- 
ample). However, in the case described the probability dis- 
tribution is Gaussian, 



g(m)dm ■ 



■. exp 



(rh-{m)f 



(7) 



where the width of the distribution grows with r. Equation 
shows that r increases with cluster temperature, which in 
turn, increases with cluster mass. Therefore, from equation 
((6]) we should expect to observe greater variability in more 
massive systems. However, the ratio {rh)/a may decrease 
with increasing cluster mass, although it is not clear how 
tcooi an d rh/rh c scale with cluster mass. 

Discretising equation ([5]) gives the s tandard form of 
a first-order auto-regr essive process (e.g. iKarlin fc Tavlorl 
ll975l ; lKing et al.ll2004 ). 



X t = 



1- ^ J." , -rr s Atf{t). 



(8) 



where At is the time-step of integration. This is a first-order 
auto-regressive process because the variable at time t de- 
pends on only the value at t — At, and is only physically 



generate Xt as a fun ction of time, variance given by (e.g. 
IKarlin fc Tavlorl 1 19751 ), 



2 

°X 



(9) 



where cj> = 1 — At/r. Note that, in this case At = 1, and 
in the limit that r > 1, a| = o\t j2, in agreement with 
equation © . Two light curves generated using equation (|17p 
are shown in figure 1. 

Even if the AGN power output is predominantly de- 
scribed by a Gaussian probability distribution, the distribu- 
tion is likely to be truncated somewhat. This is because the 
minimum AGN power is zero, while the maximum steady 
power output is governed by the Eddington limit. The trun- 
cation of the distribution will be negligible if the width 
of the probability distribution is small compared to the 
mean. However, this is probably not the case, otherwise we 
would expect to observe a high proportion of objects inject- 
ing energy at this rate. Thus, the distribution could only 
be symmetric if the Eddington limit is equal to twice the 
time-averaged AGN power output, which is approximately 
equal to Lx- A simple i nspection of the available data (e.g. 
iFuiita fe Reiprichll2004l ) demonstrates that this is probably 
not the case. 



2.1 Non-linear, asymmetric feedback 

In equation |T]) we assumed that Lx varied slowly, so that 
it could be treated as approximately constant. However, if 
the radiative losses are not matched by AGN we heating 
we would reasonably expect that the radiative losses to in- 
crease. The growth of radiative losses is a positive feed- 
back effect, so it is easy to see that the system overall 
could be non-linear and asymmetric for perturbations about 
the current, and temporar y, steady-state. For example (e.g. 
Ivon Storch fc Z wicrs 200l|), 



= Co — ciX - c 2 X 2 + f, 

at 



(10) 



meaningful if ^ 1 



^ 1. Using equation ((8)l we can 



where a > 0, and / is the white-noise random variable. 

Regardless of these considerations, given that the accre- 
tion rate is positive, definite, it seems likely that a skewed 
probability distribution will result, whatever its origin. 

In equation (|10|) . the damping is stronger to the lower 
side of the equilibrium state, than for deviations to the larger 
side. As a result the probability density function (pdf) is no 
longer Gaussian, but will be postively skewed and deviations 
to larger values are more likely than to smaller. To better 
understand this, we can think of the feedback equation as a 
randomly forced, classical particle moving in a potential, U. 
The deterministic forcing component is given by the negative 
derivative of the potential, 4j = — + /, thus for equa- 
tion ((5]): U = X 2 /(2t) is symmetric, so we would expect the 
probability distribution of Xi also to be symmetric. In con- 
trast, for equation fTO}: U = -c X + (l/2)c 1 X 2 + (l/3)c 2 X 3 . 
This is clearly asymmetric and the pdf of the Xi values will 
be skewed towards the less steep slope. 
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Figure 1. Comparison of perturbations (Xt, in model equations) generated with the same source of noise, but for different values of </>, 
where the standard deviation of the noise is unity. On the left <j> = 0.5, corresponding to Ai/r = 1/2. On the right <j> = 0.9, corresponding 
to Ai/r = 1/10. Notice that for larger <j> the range of the variability is greater and occurs on longer timescales. 



3 THE CORRELATION BETWEEN 

INSTANTANEOUS ACCRETION RATE AND 
JET POWER 

AGN jets are commonly observed to inflate cavities in the 
X-ray emitting gas that surrounds galaxies and permeates 
cluster. The jet powers can be estimated by div iding the cav- 
ity enthalpy by a characteristic timescale , (e.g. lBirzan et all 
12004 iDunn et al.ll2005l ; I Allen et al.|[2006l ). 



■Pjet — 



7 PV 



7' 



1 t 



(11) 



where P is the ambient pressure, V is the volume of the cav- 
ity, 7 is the adiabatic index of the material within the cavity, 
and t is the timescale. The pressure, and volume can be esti- 
mated from X-ray observations of galaxies and clusters. Note 
that this assumes that the cavity was inflated slowly, and in 
approximate pressure equilibrium with its surroundings. 

The appropriate timescale is not clear, and there are 
several options to choose from: the time taken to buoyantly 
rise to the observed location, the time taken to refill the 
displaced volume as the cavity rises, and the sound crossing 
time for ICM material to traverse a distance comparable 
with the bubble radius. 

The jet power, as given by equation (llip , represents a 
time-averaged quantity over the duration of the cavity in- 
flation. As a result it should corr elate well w i th th e mean 
accretion rate during this period. lAllen et al. I (|2006l ) calcu- 
late the jet powers for a sample of 9 nearby X-ray luminous 
elliptical galaxies with good optical velocity dispersion mea- 
surements. The optical velocity dispersion measurements 
permi t black hole mass es timates via the Mbh — o relation 
(e.g. Tremaine et al.ll2002l ). From this, and the X-ray data, 
I Allen et al.l ( 20061 ) also estimate the Bondi accretion rate in 
the vicinity of each supermassive black hole. It should be 
stated here that the value of the Bondi accretion rate is an 
instantaneous value which could vary significantly from the 
mean. They find the following correlation between the Bondi 
accretion power and jet power, 

log(PBondi/10 43 ergs _1 ) = 0.65 + 0.771og(P jct /10 43 ergs" 1 ) 

(12) 

Interestingly, the correlation is not linear which may sim- 
ply be a consequence of the relatively small sample, or a 
systematic error in calculating the jet power, for example. 
Alternatively, the non-linear correlation may indicate that 
other important physical effects are operating: perhaps this 
is an indication of feedback between the jet and its fuel sup- 



ply, or that the jet production e fficiency is a function of jet 
power (e.g. iNemmen et al.ll2007M . 

The possible explanation presented here is that the non- 
linear correlation arises directly from the distribution of the 
fluctuations in the instantaneous accretion rate. If the prob- 
ability distribution of the accretion rate is symmetric about 
the mean then we would expect, on average, to observe a 
linear correlation between jet-power and instantaneous ac- 
cretion rate. As a result this symmetric model of feedback 
given in equation ^ cannot provide an explanation for 
the possible non-linear correlatio n between jet pow er and 
Bondi accretion rate reported by lAllen et all (|2006T ). even 
if the variance of the distribution is a function of the mean 
accretion rate. However, a skewed probability distribution 
of accretion rates may be a possible explanation, since the 
mode and mean will not be identical. Yet, even this fact is 
not enough to explain the observed non-linear correlation: 
if P(fn/ n ^ 1) = constant 7^ 0.5, and not a function of 
H (or the variance) the correlation between jet power (/x) 
and instantaneous accretion rate m will still be linear. If 
P(rh/fj, ^ 1) > 0.5, the best-fit line determined from obser- 
vations may over-estimate the normalisation of the relation- 
ship since we expect to observe larger accretion rates more 
frequently. The reverse is also true. However, if the skew- 
ness of the probability distribution is not constant, but is 
a function of fj,, this could lead to the observed non-linear 
corr elation. 

iNipoti fc Binnevl (|2005f ) assumed that the X-ray lu- 
minosity of the cluster gas was slowly varying, while the 
AGN power varied much more rapidly. In a similar manner 
we will assume that the jet power, estimated using equa- 
tion (|lip . is also slowly varying compared to the instanta- 
ne ous Bondi accretion rat e. Following the same arguments 
as INipoti fc Binnevl (|2005l ) further, it makes sense to define 
the variable y as the ratio of the accretion power obtained 
from the instantaneous Bondi accretion rate, -Prjondi and the 
jet power P jct , 



(13) 



Again, we can argue that y is log-normally distributed, 



g{y) 



2nSy 



exp 



(lny-M) 2 



2S 2 



(14) 



for y > 0, where M and S are the mean and standard devi- 
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ation of \ny. The expectation of y is 



(y) = exp M + 



(15) 



This quantity can be directly estimated from the observed 
best-fit relation between Peondi and Pjct- Taking the Bondi 
accretion rate- Jet power relation from equation (|12|l and 
assuming that M = 0, for simplicity, we find, 



-PBondi 



-.10 n-0.23 



(y) « best - fit( 1= ) = 3.4 x 10 iu P jet "~ = , x P 



2 

(16) 

This tells us that the variance of the distribution is a func- 
tion of jet power, 



S 2 w 3 - 0.46 In I Pjct 



10 43 ergs~ 



(17) 



so that S 2 ~ 3. Such values are comparable with those found 
bv lNipoti fc Binnevl (|2005h in their model. 

The duty cycle for a source with light curve L(t) is 
defined bv lCiotti fc Ostrikerl i|200ll ) as, 



5 S ^0- = exp(S 2 ) 



(18) 



Substituting equation (|17[) into equation (|18[) tells us that 
the duty cycle is also a function of jet power (probably in- 
directly) , 



^°-° 5 (i^) ( 



(19) 



Since there is probably a close balance between the jet power 
and the radiative losses of the X-ray halo around the galaxy 
we may assume that (Pjct) ~ Lx- Therefore, we can argue 
that this result strongly suggests that AGN in larger galax- 
ies, hence larger Lx, are active for a greater fraction of the 
time. In fact, in elliptical galaxies Lx is roughly propor- 
tional to the squ are of the mass of the black hole, mbh, (e.g. 
iBest et ai1l2005l ) at the centre of the galaxy. Therefore we 
might expect S oc m^if in elliptical galaxies. 

In real it y the observed non-linear correlation given in 
[Ail en et all (2006) may not be real, but only a statistical 
anomaly. Nevertheless, it is interesting to note that the im- 
plications for the dut y cycle qualitati vely agree with the 
observational result of lBest et al.l (|2005l ) . Their result shows 
that the fraction of AGN in elliptical galaxies, which are de- 
tected in the radio, increases with the mass of the black hole 
at the centre of the galaxy. This can be interpreted proba- 
bilistically as the fraction of time the AGN is active, i.e. the 
duty cycle. At low black ho le masses (10 6 — 1 7 Mp), the ob- 
served relationship given bv lBest et all (2005) scales 
although this flattens sigificantly at higher masses. Interest- 
ingly, the estimate of the duty cycle presented here applies 
for black hole masses of ~ 8 — 9 x 10 s Mp, wh i ch agr ees qual- 
ititively with the flattening in the lBest et al.l (|2005l ) relation 
at large black hole masses. 



3.1 Stochastic variability and gentle AGN heating 

Under certain circumstances, it is possible to imagine that 
AGN activity is triggered periodically, or exhibits periodic 



activity of some sort. This could lead to large variations 
in the temperature, and thus large deviations from equi- 
librium. Yet, observations of metal distributions and cool- 
ing time profiles in clusters seem to suggest that there is 
an extremely close balance between heating and cooling in 
galaxy clusters. This is what one would expect from contin- 
uous feedback-controlled heating. However, in the simplest 
case the ensuing steady-state could result in the AGN be- 
ing permanently active with a constant power output (e.g. 
iHoeft fc Bruggenll2004l ). (This example is of interest consid- 
ering the temporal form of the accretion rate, which appears 
very similar to equation {JXJ) , in some cases). However, this 
is not what we expect from the stochastic feedback model 
discussed above. Instead, the rapid stochastic variations in 
the accretion rate mean that the AGN injects energy across 
a wide range of timescales. Thus, the AGN can maintain 
a much closer balance with its surroundings than if it was 
periodically activated. 

For a skewed probability distribution of AGN power 
output, much of this heating must occur at relatively low 
powers and so will provide the gentle heating that is con- 
sistent with observations. The power spectral density of the 
equation ([5]) described above is typical of a 'red-noise' pro- 
cess (e.g. IVaughan et al. I l2003h . with r( a 



b+ui 2 



where a 

and b are constants, and u; is the angular frequency of the 
variations. 



3.2 Stability of the system 

Time-delay in control systems tends to reduce the stability 
of the system. Indeed, for sufficiently long delays the system 
will be completely unstable. In the galaxy cluster example, if 
the heating response from the AGN is longer than the local 
cooling time, the system is unlikely to be stable. 

The main components that contribute to the time-delay 
are: 1) the time taken for inflowing gas to reach the black 
hole, and 2) the time taken for the energy fedback into the 
ICM to be dissipated. Both components are difficult to esti- 
mate, but the former moreso since it is probably composed 
of at least two stages: the inflow from large radii, and inflow 
through an accretion disc. For example, the gas presume- 
ably has som e angular mom entum, and should stagnate at 
a radius (e.g. ISarazinlll986l . and references therein), 



I 2 



r st = 



GM{< r st ) ' 



(20) 



where I is the angular momentum per unit mass of the gas, G 
is Newton's gravitational constant, and M(< r st ) is the mass 
within the stagnation radius. Th is radius is ex pected to be 
of the order of a few kiloparsecs ([Sarazi ni l 19861 ). Dynamical 
interactions could then scatter the gas clouds and allow a 
fraction to move closer to the black hole. 

The majority of the material will pass through an accre- 
tion disc that encircles the black hole; the inflow timescale 



1 In an alternative model we could replace (Lx — H) with J (Lx — 
H)dt in equation JTJ, thus the AGN power output would vary 
sinusoidally between and KiLy^i ! K\. The period of the AGN 
variability would be r = 2-k / <J\K\ K2 ) • In this case K\ would 
have a different value to the model described in the main text, 
but would still probably increase with galaxy, or cluster, mass. 
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being given by the viscous drift timescale t v i s ~ r 2 /v, where 
v is t he visco s ity an d r is the radial distance from the black 
hole. ICzernvl (|2004h gives this time as ~ lMyr at roughly 
0.1 pc for a disc around a 10 s Mq black hole. 

The time taken to dissipate the injected energy is also 
unknown, but it is concievable that the heating response 
of an AGN to a particlar event can have a significant time 
lag which could destabilise the system. This could have dra- 
matic consequences were it not for the presence of small 
stochastic fluctuations which can give rise to extremely pow- 
erful, long lived AGN outbursts. Such outbursts can occur at 
any time and may be sufficient to stabilise the system. There 
is a slight irony in this possibility since in other systems, 
stochastic variability is usually thought of as a destabilising 
force. 



4 SUMMARY 

There are many processes in the accretion of material on to 
an AGN that occur on relatively short timescales compared 
to the cooling time of the X-ray emitting gas in and around 
elliptical galaxies and galaxy clusters. In a feedback system 
these stochastic variations have been shown to lead to much 
longer and larger fluctuations in the AGN power output. 
Larger fluctuations occur in systems with longer response 
times, e.g. rich galaxy clusters. The rapid variability also 
means that power is injected across a range of timescales 
allowing an intimate balance between heating and cooling. 

The possible non-linear corr elation between B ondi ac- 
cretion rate and jet power in the I Allen et al.l (|2006l ) sample 
can be explained by a skewed probability distribution for the 
ratio of the instantaneous accretion rate to the jet power. 
This explanation also predicts that the duty-cycle of AGN 
activity varies with its environment, or black hole mass. 

Stochastic variability of the accretion rate may help to 
stabilise the system against the large time-delays that are 
probably inherent with AGN feedback. 
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